Flat bands with local Berry curvature in multilayer graphene 
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We demonstrate that flat bands with local Berry curvature arise naturally in chiral (ABC) mul- 
tilayer graphene placed on a boron nitride (BN) substrate. The degree of flatness can be tuned by 
varying the number of graphene layers N . For N = 7 the bands become nearly flat, with a small 
bandwidth ~ 3.6 meV. The two nearly flat bands coming from the K and K' valleys cross along lines 
in the reduced zone. Weak intervalley tunneling turns the bandcrossing into an avoided crossing, 
producing two nearly flat bands with global Chern number zero, but with local Berry curvature. 
The flatness of the bands suggests that many body effects will dominate the physics, while the local 
Berry curvature of the bands endows the system with a nontrivial quantum geometry. The quantum 
geometry effects manifest themselves through the quantum distance (Fubini-Study) metric, rather 
than the more conventional Chern number. Multilayer graphene on BN thus provides a platform 
for investigating the effect of interactions in a system with a non-trivial quantum distance metric, 
without the complication of non-zero Chern numbers. We note in passing that flat bands with 
non-zero Chern number can also be realized by making use of magnetic adatoms, and explicitly 
breaking time reversal symmetry. 



The behavior of strongly interacting electrons and the 
effect of quantum geometry are two of the most ex- 
citing fields of research in modern condensed matter 
physics. Strong interactions can give rise to effects such 
as fractionalization [iHS], where the elementary excita- 
tions carry some fraction of the 'fundamental' electronic 
quantum numbers. One popular way to access strong cor- 
relation physics is to consider band structures that are 
flat (non-dispersive) , since in such band structures inter- 
actions naturally dominate over kinetic energy. Mean- 
while, the two ingredients of quantum geometry are in- 
tegrated Berry curvature (a.k.a. Chern number for two 
dimensional systems) and Fubini-Study metric U, the 
metric measuring the quantum distance. While non-zero 
integrated Berry curvature can lead to the existence of 
states that are protected against disorder [SHT] and to ef- 
fects like quantized Hall conductance of a band insulator 
[5], non-zero Fubini Study metric can lead to phenom- 
ena such as pseudospin conservation laws in single layer 
graphene j4j and unusual features in the current noise 
spectrum of a band insulator |9]. Recently, a strong in- 
terest has emerged in flat band systems as playgrounds 
for investigating the interplay of strong correlation effects 
and non-trivial quantum geometry. However, attention 
has been mostly focused on systems with non-zero inte- 
grated Berry curvature, [TTI - 120] . while the effect of the 
Fubini-Study metric has largely been ignored. 



The influence of quantum geometry on the strong cor- 
relation physics stems from the fact that the interactions 
couple to electron density, and the electron density op- 
erators projected onto the flat band obey a non-trivial 
commutation relation if the band has a non-trivial quan- 
tum geometry [IBIET]. The general commutation relation 
for projected densities in a band with non-uniform Berry 



curvature is [5T]: 
[Pqi,Pq2] ~ iqiAq2^[S(k)u;;(k+)ufa(k_)7^_^7k_] (1) 

where ql A q2 = z.(qi x qa), k± = k ± ai±a2 ^ ^j-^) 
is the local Berry curvature for the band of interest for 
the single-particle hamiltonian H = a b '^k a^a.bO^)cii,b 
and the corresponding eigenstate is given by |k) = 

7k |0) ~ Sa"a(^)'^ka l'^)' Thus far, attention has been 
focused on systems where the Berry curvature is nearly 
uniform pTk20| . However, it is apparent from ([T]) that 
non-trivial quantum geometry effects do not require a 
uniform Berry curvature. The influence of the quantum 
distance metric may be most clearly revealed in a flat 
band system with local Berry curvature, but zero inte- 
grated Berry curvature, since here the non-trivial quan- 
tum geometry effects (encoded e.g. in the projected den- 
sity commutator) arise purely due to the metric. Thus 
far, theoretical studies have largely ignored this excit- 
ing direction of research, in part because of the lack of 
experimental realizations. 

In this Letter, we show that flat bands with local 
Berry curvature (but vanishing integrated Berry curva- 
ture) arise naturally in chiral multilayer graphene. Our 
proposal exploits the fact that ABC stacked multilayer 
graphene in the presence of a perpendicular electric field 
has a bandstructure with flat pockets that possess Berry 
curvature. Placing the graphene on a hexagonal Boron 
Nitride (BN) substrate then produces a superlattice po- 
tential such that the reduced Brillouin zone lies entirely 
within the flat pocket. Umklapp scattering opens a 
bandgap at the reduced zone edge. For a TV layer sys- 
tem with N > 5 layers, the lowest band is nearly flat, 
with a bandwidth 5meV. We have verified that this 
nearly flat band has a non-vanishing local Berry curva- 
ture. Chiral multilayer graphene thus provides an ideal 
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platform for investigating the interplay of strong correla- 
tions and quantum geometry, with the quantum geome- 
try effects coming from the (hitherto neglected) channel 
of the Fubuni-Study metric, rather than the more con- 
ventional channel of non-zero Chern number. As such, 
we believe chiral multilayer graphene offers an exciting 
new frontier in the study of systems with an interplay of 
strong correlations and quantum geometry. 

We note in passing that a non-zero integrated Berry 
curvature can also be obtained by using magnetic 
adatoms rather than vertical electric field to open a gap 
between conduction and valence bands. 

ABC stacked graphene: A single graphene sheet con- 
sists of a honeycomb lattice of carbon atoms. The 
honeycomb lattice consists of two sublattices, A and 
B. Chiral multilayer graphene consists of graphene 
sheets with an ABC stacking order (each succeed- 
ing sheet is rotated by 27r/3 relative to the preceding 
sheet). Every lattice site in the bulk is either directly 
above or directly below another lattice site. In the 
('0iAkV'iBkV'2yik'02Bk i^NAkipNBk) basis, the near- 
est neighbor tight binding Hamiltonian for an N layer 
system takes the form [22 
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Here, tp = io(exp(ifc^a) 2exp(-i^) cos(^%i^)), 
represents nearest neighbor hopping within each 
graphene layer {to ~ 3eV), and 7 ~ 300meV represents 
interlayer hopping between two sites that lie on top of 
each other. Here a is the lattice constant of graphene. 

In the absence of interlayer hopping, 7 = 0, the band- 
structure consists of A'' copies of the graphene bandstruc- 
ture, E — ±\tp\ [IS]. Near the two inequivalent corners 
of the Brillouin zone, conventionally labelled K and K' , 
the function tp vanishes as tp w wp_|_ and tp w vp 
respectively, where p± = p^, ± ipy and v — 3ato/2 is the 
Fermi velocity for graphene. 

We now consider interlayer hopping 7 7^ 0. This inter- 
layer hopping causes all the bulk sites to dimerize, open- 
ing up a bulk gap of size 7 at the Dirac points. On the 
top and bottom surfaces, there are undimerized lattice 
sites, which give rise to gapless surface states. The low 
energy single particle Hamiltonian for the surface states 
of an A^ layer graphene takes the form [21 [Ml UHl 123 
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Here Hk is the Hamiltonian in the K valley, and 
the Hamiltonian in the K' valley is given by Hk'{v) = 
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FIG. 1: Low energy band structure for A'^ layer ABC stacked 
graphene in the presence of a vertical electric field. The band 
structure is plotted in the vicinity of the K point, assum- 
ing that the potential difference between the top and bottom 
layers A = 0.167fo ~ SOmeF. 



i?J^(— p). The basis is such that (1, 0) is a Bloch state in 
the A sublattice of the top layer, and (0, 1) is a Bloch 
state in the B sublattice of the bottom layer. Only 
the lowest order terms have been written in each ma- 
trix element. This Hamiltonian is valid up to an energy 
scale 7. Note that this energy scale is independent of 
A^, the number of layers of graphene. The above effec- 
tive Hamiltonian consists of a single bandcrossing with 
non-trivial Berry phase A^tt, and low energy dispersion 
E = zLv^p'^ /'^^~^ . Thus, the conduction band has a 
very flat pocket of a size p ~ 7/1; around the K and K' 
points, and this pocket becomes perfectly flat in the limit 
A^ — >■ 00, corresponding to 'rhombohedral graphite'. The 
emergence of this flat pocket has previously been dis- 
cussed in [Hm]. 

Now consider the application of a vertical electric field 
leading to potential difference V between any 2 consec- 
utive layers. We assume that 7 :s> e{N — 1)V. This 
ensures that the surface bands are much below the bulk 
bands. The low energy effective Hamiltonian then takes 
the form 




(4) 



The resulting band structure is plotted in Fig[T] We can 
see that a gap of magnitude A = (A^ — l)eV opens in 
the surface state spectrum. We can also see that the 
band near the K point is extremely flat, asymptoting to 
perfect flatness in the limit N 00. 

Effect of BN substrate: BN is a popular substrate for 
graphene. It has the same hexagonal structure, but with 
a slight lattice- mismatch (lattice constant ~ ||a) and is 
a large band-gap insulator {Egap ^ 5eF). Placing multi- 
layer graphene on BN substrate introduces a superlattice 
potential with a. Z x Z supercell, where Z — 56. The 
primitive lattice vectors corresponding to the reciprocal 
lattice are 
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FIG. 2: Bandwidth A of the lowest conduction band for 
layer chiral graphene. For N > 5, the bandwidth comes 
mainly from umklapp scattering at the zone boundary. 



The strong superlattice potential is seen primarily by 
states on the lower surface of the multilayer graphene, 
close to the substrate. These states are either conduc- 
tion band or valence band states, depending on the sign 
of the vertical electric field. We assume the electric field 
is chosen such that it is the conduction band that mainly 
sees the superlattice potential. If the multilayer graphene 
were sandwiched between BN sheets, then both conduc- 
tion and valence bands would see a strong superlattice 
potential, irrespective of the sign of electric field. 

The reduced 'Brillouin zone' for the superlattice is 
hexagonal, but scaled by a factor 1/Z with respect to 
the original Brillouin zone. If Z is chosen such that 
\K\/Z < j/v, i.e. Z > Zc, where 



Zc = 207r/\/3 = 36.3 



(6) 



then the 'flat pockets' at K and K' extend over the entire 
reduced Brillouin zone. Moreover, umklapp scattering at 
the boundaries of the reduced Brillouin zone opens a gap 
between the flat pocket and the rest of the band. Since 
for BN substrate, Z = 56 > Zc, it follows that placing 
chiral multilayer graphene on BN leads to flat bands that 
extend over the entire reduced Brillouin zone, and which 
are separated from the rest of the bandstructure by the 
umklapp energy scale A. 

We now explicitly calculate the bandstructure in the 
reduced zone. It may be readily determined that for the 
above superlattice potential, the K and K' points in the 
original zone get folded to inequivalent corners K and K' 
of the reduced zone. Meanwhile, the reciprocal lattice 
vectors of the superlattice satisfy IB1.2I < 7/w < 2|Bi.2|, 
where 7/u is the width of the flat pocket. Since Bragg 
scattering is only effective between states that are near 
degenerate, we restrict our attention to Bragg scattering 
by a single reciprocal lattice vector, and neglect higher or- 
der Bragg scattering events. The matrix elements of the 
superlattice potential (k|y |k + Bi^2) were determined 
by modelling the BN superlattice as a positive (5-function 
potential at each B site and a negative 5-function poten- 
tial at each N site [H] . Hence we obtained the bandstruc- 
ture in the reduced zone. Although the bandstructure 



contains Z'^ = 56^ = 3136 conduction bands per spin and 
valley, only the lowest conduction band is flat, and is sep- 
arated from the higher bands by an energy scale of order 
A, where the gap may be estimated from the DFT cal- 
culations in [57], and is of order lOmeV. Meanwhile, the 
bandgap A between conduction and valence bands can 
be externally controlled using gates, and may be made 
as large as desired. We should take A :s> lOmeV to en- 
sure that we do not mix conduction and valence bands. 
For specificity, we suggest using A — 50meV, which is 
easily achievable by gating [55]. 

The superlattice potential from the BN substrate also 
introduces intervalley tunneling of magnitude X/Z [29j . 
which turns the bandcrossings of the two nearly flat 
bands coming from K and K' valleys into avoided cross- 
ings. The resulting bandstructure contains two strictly 
non-degenerate flat conduction bands, separated by a 
minigap X/Z 0.2meV. The lowest energy band comes 
mostly from the K valley in regions closest to the K cor- 
ner of the reduced zone, and comes mostly from the K' 
valley in regions closest to the K' corner. From the band- 
structure calculations, we can extract the bandwidth of 
the two low-lying conduction bands. In Figj2j we plot the 
bandwidth of the lowest conduction band as a function of 
N. We see that for TV > 5 the two low-lying conduction 
bands are nearly flat, with a small residual bandwidth 
around 5 meV which comes mainly from umklapp scat- 
tering at the zone boundary. The N = 7 layer system 
actually has minimum bandwidth of ^ 3.6 meV, due to 
a cancellation between 'intrinsic' curvature of the band 
and the effect of umklapp scattering. 

We also note that there is an indirect band gap sepa- 
rating the two flat bands from the higher energy non-flat 
bands. It vanishes for TV < 6 and then decreases with 
increasing N. Transitions across the indirect band-gap 
must be phonon assisted and should be weak, but never- 
theless it is essential that we work with > 5 to have a 
truly isolated fiat band. Fortunately, A = 7 is optimal 
not only through having the fiattest band (Figj2]), but 
also because its indirect band gap is of order the direct 
band gap. For N — 7 the direct band gap is ^ 6meV. 

We have now verified that the lowest conduction band 
can be made flat. Now we show that this band has non- 
vanishing Berry curvature. The momentum-space Berry 
curvature for the nth band is given by [8 

B (k) - ^ 2/to {nk\v.^\n'k) (n'k|i;^|nk) 



where is the velocity operator and £"„ is the eigenen- 
ergy corresponding to the |nk) eigenstate. It may be 
readily determined that the K and K' bands have oppo- 
site signs of local Berry curvature. The conduction and 
valence bands also have opposite signs of local Berry cur- 
vature at every point in momentum space. As a result, 
when the chemical potential is placed between conduc- 
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FIG. 3: Contour plot of Berry curvature in the lowest flat 
conduction band for N = 7. The red/yellow regions come 
mainly from the K valley and have positive curvature, while 
the blue regions come mainly from the K' valley and have 
negative curvature. The Berry curvature integrated over the 
band is zero. 



tion and valence bands, the system displays quantum 
valley Hall effect [SU]. However, flat band physics will 
manifest itself when the chemical potential is placed in- 
side either the conduction or the valence bands. In this 
case we can focus on the band that contains the chemical 
potential. 

We consider the conduction band for specificity. Due 
to avoided crossings between bands coming from the two 
valleys, the lowest conduction band contains regions with 
positive and negative local Berry curvature respectively 
while the integrated Berry curvature (Chern number) is 
zero (Figjsj). Meanwhile, there is a second nearly flat con- 
duction band which is separated by a mini gap equal to 
the inter valley scattering amplitude. Interactions may 
enhance this minigap, lifting the near degeneracy of the 
two conduction bands to gain exchange energy, in accor- 
dance with 'quantum Hall ferromagnetism' [lOj . 

Thus, we conclude that placing N — 7 layer ABC 
graphene on BN substrate allows us to realize two nearly 
flat bands with local Berry curvature but zero Chern 
number which are separated from the non-flat bands by 
a band gap of GmeV and are separated from each other 
by a minigap of 0.2meV, equal to the inter valley scat- 
tering amplitude. Chiral multilayer graphene thus of- 
fers a promising playground for investigating the effect of 
strong interactions in the presence of a non-trivial quan- 
tum distance metric. 

Chern number from adatoms. The key feature of chi- 
ral multilayer graphene is that it allows access to a sys- 
tem with a non-trivial quantum geometry but without 
Chern number. However, a non-zero Chern number may 
also be obtained by making use of adatom deposition on 
the outermost graphene layers to open up a gap between 
conduction and valence bands, rather than using vertical 
electric field. Adatom deposition also introduces a su- 
perlattice potential. For the appropriate choice of (time 
reversal symmetry breaking) adatoms, the two valleys ac- 
quire the same sign of Berry curvature |32H34| . The rest 



of the analysis proceeds exactly as before, only now the 
Berry curvature has the same sign everywhere in the flat 
band, and thus does not cancel. 

Concluding remarks. We have shown that ABC 
stacked multilayer graphene placed on BN substrate has 
a bandstructure containing flat bands. The seven layer 
system is ideal for this purpose. The flat bands have 
nonzero local Berry curvature but zero Chern number. 
Thus, chiral multilayer graphene represents an exciting 
new frontier in the study of interaction effects in sys- 
tems with non-trivial quantum geometry, allowing access 
to an interaction dominated system with a non-trivial 
quantum distance metric but without the complication 
of a non-zero Chern number. 
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SUPPLEMENT 



A. We consider the case of multilayer graphene with interlayer hopping 7^0. It is con- 
venient to separate the low energy and high energy subspaces by writing the Hamiltonian in the 
('0iAk'/'A'SkV'isk'*/'2yik'02Bk'03Ak ''P{N~i)Bk4'NA'k) basis. The Schrodinger equation in this basis becomes 



H21 H22) yiphighj \^high 



(8) 



Here e is the energy eigenvalue, ipiow is a two component spinor, iphigh is a 2A'' — 2 component spinor, Hn is a 2x2 
zero matrix and H12 and H22 are 2 x (2N-2) and (2N-2) x (2N-2) matrices respectively. Here H12 and H22 have the 
form below 
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We can eliminate the high energy dimer components in the system of equations ([8]) to get a Schrodinger equation for 
only the low energy components in the manner of ^26]. Since we are interested in only the low energy states, we make 
a binomial expansion to the first order in e assuming it is very small compared to H22- We get (-ffe// ~ £)4'iow = 
where H^ff = {I + Hi2{H22)~^ H2i)~^ {Hu — Hi2{H22)~^ H21) ■ This leads to a low energy single particle Hamiltonian 
for the surface states of an N layer graphene which takes the form 
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B. First, we calculate the matrix elements of the external potential V coming from the BN between Bloch states, 
(kjFjfc + q). This is just equivalent to calculating the Fourier transform of the BN potential. We model the BN 
potential as X'5{v{) — A'i5(rj), where ri are the positions of boron atoms and rj are the positions of nitrogen atoms. 
This is a natural model, since the atomic numbers of B and N are one less than and one more than carbon respectively. 
Taking the delta functions to have slightly different weights will not qualitatively alter our results. The boron atoms 
sit on the A sublattice of the hexagonal superlattice, and the nitrogen atoms sit on the B sublattice. 

The Fourier transform of the above potential takes the form 



(k|y|k + q) =l/(q) = 5k,Q/i(q) 



(12) 



where Q denotes a reciprocal lattice vector of the BN lattice (which is equal to Bi_2 in Eq|5]of the main text modulo 
reciprocal lattice vectors of the graphene lattice) and /i(q) is a form factor coming from the two site nature of the BN 



unit cell. For the model potential under consideration, /i(q) oc (e''^ '' 
(k|l/|k + B2) = i\, (k|l/|k + Bi B2) = zA, and /i(-k) = A*(k). 



e^qi-j). Hence, we obtain (k|y|k-|-Bi) 



6 



Thus far we have assumed that the BN superlattice potential can be modelled as a delta function array. In fact, 
the B and N atoms carry + and — charge respectively. The BN superlattice potential may thus be better modelled as 
a delta function array convolved with a 1/r envelope. The 1 jr envelope simply reflects the Coulomb potential arising 
from a local charge imbalance. This Fourier transforms to a delta function array multiplied by a 1/fc envelope. Thus 
intervalley scattering by Z reciprocal lattice vectors is weaker than intravalley scattering by one reciprocal lattice 
vector by a factor of \ jZ. However as we will see below, this weak intervalley scattering is still important along lines 
in the reduced zone where the K and K' bands are degenerate. 

Now we calculate the band structure. First we consider the low energy hamiltonian of multilayer graphene without 
BN substrate in Eqj4j It is written in a basis such that (1, 0) is a Bloch state in the A sublattice of the top layer 
and (0, 1) is a Bloch state in the B sublattice of the bottom layer. Now we consider the case of multilayer graphene 
sandwiched between BN sheets. Without considering the folding of the bands at the reduced zone edges, the 4 low 
energy bands corresponding to a particular hexagonal unit cell with center Q are found from the eigenvalues of 
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where Q' points to the K point of the hexagon, for which the distance between k and Q' is the minimum. A similar 
definition holds for Q" and K' point. The matrix is written in a basis where (1, 0, 0, 0) denotes a state on sub lattice 
A and with valley (0,1,0,0) denotes a state on sub lattice B and valley (0, 0, 1, 0) denotes a state on sub lattice 
A and valley K' and (0, 0, 0, 1) denotes a state on sub lattice B and valley K' . In the off diagonal blocks wc have 
included a weak inter-valley scattering coming from the superlattice potential. 

The band folding at the zone edges can be taken into account by considering the four low energy bands arising 
from a central hexagonal unit cell in the reciprocal space and those arising from it's 6 nearest neighbouring cells. The 
corresponding band structure is easily found from the eigenvalues of the following 7x7 matrix: 
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where Kq points to the center of the central hexagon, Ki, K2, K3, K4, K5 and Ke point to the centers of the 6 
adjoining hexagons. The matrix is written in a basis where (1, 0, 0, 0, 0, 0, 0) is a state near , (0, 1, 0, 0, 0, 0, 0) is a state 
near K\, (0, 0, 1, 0, 0, 0, 0) is a state near K2 etc. and the matrix corresponding to the intra- valley scattering between 
the valleys of neighbouring hexagonal unit cells is 
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This 4x4 matrix is written in the same basis as Eq.l3, where (1, 0, 0, 0) denotes a state on sub lattice A and with 
valley K, (0,1,0,0) denotes a state on sub lattice B and valley K, (0, 0, 1, 0) denotes a state on sub lattice A and valley 
K' and (0,0,0, 1) denotes a state on sub lattice B and valley K' . 



